Nonlinear wave interaction and spin models in the MHD regime 



G. Brodin, J. Lundin, J. Zamanian and M. Stefan 
Deptartment of Physics, Umea University SE 901 87 Umea, Sweden 

Here we consider the influence on the electron spin in the MHD regime. Recently developed 
models which include spin-velocity correlations are taken as a starting point. A theoretical argument 
is presented, suggesting that in the MHD regime a single fluid electron model with spin correlations 
is equivalent to a model with spin-up and spin-down electrons constituting different fluids, but where 
the spin-velocity correlations are omitted. Three wave interaction of 2 shear Alfven waves and a 
compressional Alfven wave is then taken as a model problem to evaluate the asserted equivalence. 
The theoretical argument turns out to be supported, as the predictions of the two models agree 
£SJ , completely. Furthermore, the three wave coupling coefficients obey the Manley-Rowe relations, 

which give further support to the soundness of the models and the validity of the assumptions made 
in the derivation. Finally we point out that the proposed two-fluid model can be incorporated in 
standard Particle-In-Cell schemes with only minor modifications. 



• PACS numbers: 52.25.Xz, 52.35.Bj, 52.35.Mw 



I. INTRODUCTION 



Considerable interest has recently been devoted to the study of quantum plasmas^see e.g. Refs [1 jofl. Much of the 
research has been motivated by applications to e.g. quantum wells 0], plasmonics [ic) . spintronics astrophysics 
[l2| and ultra-cold plasmas Two of the most basic and much studied quantum effects are those of the Fermi 

pressure and the particle dispersive effects (directly associated with the Bohm De Broglie potential), see e.g. Refs. 
Q-T HI 01 • Other studies 0, [T3 - [2l| focus on the electron spin properties that result in the magnetic dipole force and a 
r/} , magnetization current, in addition to some more complex aspects of the spin dynamics. Although most quantum 
O ■ effects has a tendency to be more important in plasmas of high density and low temperature, the regimes of relevance 
" c/3 ' differ to some extent for the various quantum effects, see Ref. [l9| for a discussion of this issue. A consequence is 
that it is possible to focus on certain of the quantum effects and ignore the others. In this paper we will make use 
' of this fact and concentrate on the physics associated with the spin coupling in the Pauli Hamiltonian. Our starting 
i point is a recenty presented spin- fluid model (l6| . derived from kinetic theory [l5| . which in addition to the most basic 

spin precession dynamics includes effects of spin-velocity correlations. Evaluating this model in the MHD regime, we 
make a conjecture based on certain teoretical arguments; That there are two equivalent ways to model spin-MHD 
^ ■ dynamics, either by a one-fluid model including spin-velocity correlations, or by a two-fluid model without spin- 
velocity correlations. In the latter case the spin-up and spin-down states relative to the magnetic field are regarded 
as different fluids The conjectured equivalence of these models in the MHD regime is tested by considering a 

' specific problem of three-wave interaction. For this purpose we calculate the coupling coefficients between two shear 
• , Alfven waves and one compressional Alfven wave (fast magnetosonic wave) in a magnetized plasma. The coupling 
£!Q ■ coefficients turns indeed out to be identical in the two cases, and the coefficients are also seen to obey the Manley-Rowe 
relations, which give further support to the soundness of the models used. The applicability of the two-fluid model 
I ■ without spin-velocity correlations in the MHD -regime, which is strongly supported by our findings, is a very useful 
• • \ result. The reason is that as this model can be easily adopted into standard Particle-In-Cell schemes with only small 
. modifications, as will be discussed in the final section. 

X: 

d ■ II. MODEL EQUATIONS 

Starting from a scalar kinetic equation for a spin-I/2 particle [15) . spin fluid equations can be derived [l6| . These 
are given by the continuity equation 

cW s) + V • {n^v^) = (1) 

and the fluid momentum equation 

where the superscript s = e,i denotes the species (electrons or ions), D/Dt = dt + v( s ) ■ V, and summation over 
repeated indices i, j, k = x,y,z is implied. Here is the mass, is the charge, fj,^ the magnetic dipole moment, 
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n( s ) is the number density, and v( s ) is the fluid velocity of species s. Furthermore, S is the spin vector normalized 
to unity, is the pressure tensor, and is the Levi-Civita symbol. Since the ions normally have a much smaller 
magnetic moment than the electrons [2^, the spin contribution due to the ions can be neglected compared to the 
electron contribution, i.e. we may let fjfi> re 0. We have here neglected contributions to the force from the Fermi 
pressure and particle dispersive effects (the so called Bohm de Broglie potential), see @. For a discussion of the 
parameter regime of importance of various quantum effects, see e.g. Refs. fl9H23j. In this approximation, the pressure 
moment satisfies the evolution equation 

Dt ~ lk dx k * k dx k ij dx k + niW mn * m n + mW jmn im n 
^ 8B k M « v dB k 

where again the last two terms can be dropped for the ion species. 

Furthermore, to describe the spin dynamics we need the electron spin evolution equation, which is given by 

Dt ~ h 3 3 Uk m (e) n (e) ft^ W 

where is the spin-velocity correlation tensor. Finally the evolution of the spin- velocity moment is described by 



~1h£ e) , B , a d 



s (e) uv k _ ^—1 P (e) - 



Dt ij dx k ik dx k jk dx k 

\ £ jkl^ik a l H £ — e ikl^ki °l + M n q M n b i h k "5 W 



In Eq. © we have neglected the heat flux tensor Qijk to obtain a closed set of equations. Similarly we have neglected 
the higher order tensor Aijk in the evolution equation for the spin- velocity tensor Eq. ((5|). The validity of the 
truncation has been investigated in Refs. (l5l. l24l|. and the truncation seem to be an accurate approximation in the 
low-temperature limit. The equations above together with Maxwell's equations constitute a closed system, where a 
magnetization current density ]m — V x M, due to the spin, should be added to the free current density, and where 
naturally all species contribute in the latter term. The set of equations (HI)-© has been studied by Refs. [THHH; but 
without inclusion of the ion dynamics. The aim of the current paper is to apply the above set of equations to the 
MHD regime where the ion dynamics is essential, at the same time as making a careful evaluation of the electron spin 
magnetization. 

A. MHD-limit 

As concluded in the previous section, we will primarily consider wave dynamics in the MHD regime where the 
frequencies are smaller than the ion-cyclotron frequency and the wavelengths are longer than the Larmor radius. 
Furthermore, we will consider the low-temperature (i.e. low-beta) limit, where the pressure terms are dropped. 
Under these assumptions, the system will be described by the magnetohydrodynamic equation (20| 

p ( + u • v) u = -V (^— - M ■ + (B • V)M - VP e (6) 



together with the equations 



and 



dt J \2fJ-o 

d t B = V x (u x B) . (7) 



^+V-(pu) = (8) 



Here we have neglected the electron contribution to the fluid mass density by setting p re m^ny', and the fluid 
velocity can be written as u = (n^m^v^ + nWm'''v''))/ p re v". We have neglected the magnetic moment of 



3 



the ions such that the magnetization M = fi^n^S^ is purely due to the electron spin. The derivation of Eq. © 
made in Ref. [2(| was done starting from a somewhat less elaborate set of equations, not including the spin-velocity 
correlations. However, the derivation does not involve the Eqs. ((J)-© describing the spin dynamics, and hence we 
may adopt this result within the current model. 

Without the magnetization M we obtain standard ideal MHD equations, and thus ©, (0) and (J5| constitute a 
closed system. The magnetization can be easily included with the spin determined by Eqs. (j4| together with ([5]), where 
the terms containing derivatives of the velocity turns out to be negligible. The aim is to solve for the magnetization 
in terms of the magnetic field, in which case Eqs. (J5]), ([7]) and © are sufficient to produce a closed spin-MHD theory. 
This can be achieved in two different ways; either by considering the electrons in spin up and spin down states 
relative the magnetic field as two separate fluids, or by treating them as a single fluid with a macroscopic spin that 
is proportional to the difference in population density of the two spin states. We will now go on to discuss this in 
further detail. 



B. One-fluid vs two-fluid 

In this sub-section we will discuss how to determine the magnetization, in order to use Eqs. (jSJ), ([7]) and ([5]). To find 
the magnetization we first need to solve © to determine the spin. The first term of the right hand side of Eq. ([J} is 
the basic spin precession. If the spin-velocity correlations in @ can be neglected, the solutions for S are particularly 
simple in the MHD regime. This is because the left hand side term of Eq. ^ is smaller than the spin precession 
term by a factor of the order 0(a/ ch ) /uiig), where u/ ch ) ~ <9< is a characteristic frequency scale of the problem, and 
^cg 1 ^ ~ 2/i e B/H is the characteristic spin precession frequency (which is close to the characteristic cyclotron frequency 
w' * 1 ' ~ qB/m). Assuming that spin- velocity correlations can be omitted, we note that spin evolution equation in the 
MHD regime reduces to 

e^ e) B k = 0. (9) 

This has two solutions, where S is either parallel or antiparallel to B , that is Si = ±&j, where bi = B{j B is a unit 
vector in the direction of B. A comparatively general way to deal with spins obeying Si = ±6, is to consider a 
two-fluid model of electrons, where for one of the species the electron spin state is parallel to B, and for the other 
species antiparallel. Eq. © then implies that these spin states are conserved. However, as seen from the above 
discussion this is only an adequate approximation if the spin-velocity correlations give a small contribution in Eq. 
((3]). Thus our next step is to outline the solutions of Eq. using MHD approximations, in order to determine the 
contribution from £Vv in (QJ. Firstly we note that the three first terms in ^ are at most of order u/^Ey, whereas 
the fifth and sixth are of order ^ ch ' Ey ~ ^>cg^ Ey . Thus neglecting the three first terms we can formally write Eq. 

(HI) on the form O ■ , where O is a 9 x 9-matrix where all coefficients are ±w CQ , or ±u) cga . Here E is a 

9-component vector containing all elements of £y, and ^ is a 9-component vector containing the source terms, i.e. 
Pjk(dSi/dxk), fin{dBi/dxj) and /inSiSk(dBk/dxj), uj ca = qB a /m and uj cga = 2[i e B a /h with a = x,y,z. Note that 
the full field strength is used and not the linearized field when defining uj ca and u) cqa - Since O contains no operators, 

^ ^ y 

we can do a simple matrix inversion to find £ = O ~ l a . This turns out to be sufficient to determine all components 
of £jj, except for a component directed as b ® b. Thus this approximation scheme allows us to compute Ey apart 

from a contribution that can be expressed as Ey = Qbibj, where $ is a scalar field. Using the solution 

we can easily check that the determined components of E,y are of order Ey ~ (isnkB/tJc ■ This means that the 
contributions from Ey are sufficiently small to be neglected in (U). However, in general we must also account for 
the contribution Ey = whose magnitude is unknown. Since the scalar field 4> cannot be determined if the 

three first terms of ([3]) are omitted, we must extend our model to the solve the full case of Eq. ([5]). We will do so 
within a one-fluid model in the section III A, and it turns out that <J> becomes sufficiently large for this component 
of Ejj to signficantly influence the solutions to also within the MHD-regime. However, it also turns out that in 
order to get a large value of $, we must have a spin vector that is different from ±b. This is the normal case in a 
one-fluid theory, where the macroscopic spin results from averaging over all spin states. However, within a two-fluid 
MHD model (treating spin-up and down states as different species) without spin-velocity correlations Eq. (0) can be 
applied leading to Si = ±&i, and the situation would then again be modified. Indeed, contracting Eq. ([5]) with bibj to 
compute the source terms for $, we find that all the source terms for Ey vanishes if Si = ±6j, as the fourth term in 
Eq. ([5]) becomes bibjPjk{dSi/dxj), which is zero as bi(dSifdxj) = ±(l/2)d(bibi) /dxj , whereas terms 7 and 8together 
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become 



dBi dB k 
bibj-^ bibjSiSk-^. — = 

OXj OXj 

where Si — ±6i was used in the last step. This result provides the theoretical basis for adopting a two-fluid model 
of electrons in the MHD- regime and omitting spin- velocity correlations in (TJ|, leading to Si = ±6j. The division into 
two fluids leaves the rest of the basic equations structurally unaffected, but we now obatin two contributions such that 
the magnetization is calculated as M = /m-|-s-|- + pn^s^ due to the difference in density perturbations of the two spin 
states. We will consider this in more detail within perturbation theory in our model problem below. The conclusions 
of this section is then confirmed, since the one-fluid models that keeps the spin-velocity correlations in Eq. (U) give 
indeed an identical expression for the magnetization as the two-fluid model with up and down spins Si = ±6». The 
allowance for independent density variations of the two species in the latter model provides the physical mechanism 
that reproduces the effects of spin-velocity correlations in the one-fluid model. It should however be stressed that this 
conclusion is limited to the MHD regime. 



III. THREE WAVE INTERACTION - A MODEL PROBLEM 



We will now consider a model problem with the purpose of testing our conclusions about the similarities between 
the one-fluid and two-fluid models outlined in the previous section. Specifically, we consider three wave interaction 
between two shear Alfven waves (A, A') and one compressional Alfven wave (MS); MS — > A + A' . Using three-wave 
interaction as a model problem has the advantage that an unphysical assumption (or an incorrect calculation) is likely 
to result in a broken Manley-Rowe symmetry |25| . in which case one gets a clear indication that something needs to 
be revised. 

The waves are assumed to be small perturbations on a homogeneous background, and we write B = B$i + Bi, 
p = po + pi, etc., but omit the index I on variables whose backgound values are zero. Furthermore, we omit index 
1 whenever the cartesian components are specified for notational convenience, i.e. we write Bi = B x x + B y y + B z z. 
Moreover, we assume that there is no drift so that Uo = and also that there is no spin- velocity correlation in the 
background distribution, i.e. T,q = 0. For simplicity we further assume that the temperature is sufficiently low such 
that the equilibrium pressure can be neglected, Po — (i.e. that we have a low-beta plasma with the ion-acoustic 
velocity much smaller than the Alfven velocity). Furthermore, assuming that pbBq/ '(ksT) <C 1 we can make the 
approximation that the equilibrium spin up and spin down populations are equal so that nof = TI04. i n the two-fluid 
model which implies that the total zeroth order magnetization vanishes (2(|. For consistency between the one-fluid and 
two-fluid models we should consequently pick Mo = also in the latter case. The difference in the model equations 
between the one-fluid and two-fluid approach is then primarily that in the two-fluid model we have So = ±z for the 
two spin states (in which case we obtain a finite zero order magnetization if only one of the electron fluids are counted) 
whereas in the one-fluid model Mo = and So = 0. 

Next we make a harmonic decomposition <9 f — > —icu and <9 X — > ik for each wave, where the frequencies and wave 
vectors satisfy the conditions 

u MS =u A + u A ' (10) 

k MS = k A + k A' (n) 

with the index MS denoting the compressional Alfven (or fast magnetosonic) wave, and the index A and A 1 denoting 
the shear Alfven waves. The coordinate system is defined so that the z-direction points in the direction of the 
unperturbed magnetic field, B = Bqz, and for simplicty we assume all wave vectors to lie in the xz-plane. 

Throughout the calculation we will use uj/lj c and fcC^/wc as small expansion parameters (where uj c = qBo/m is 
the cyclotron frequency), in accordance with standard MHD theory. Here Ca = {Bq / ' pomtno) 1 ^ 2 = (Bq/ papo) 1 ^ 2 is 
the Alfven velocity, and u> and k represents any of the wave frequencies or wave vector components. We also note 
that uj cg ~ u c , where uj cg = 2p e B /h is the spin precession frequency. Furthermore, of the same order as 

the ion-cyclotron frequency, and is therefore much larger than wave frequencies within the MHD regime (27| . We will 
therefore drop terms proportional to (to C g — ^c)" 1 compared to lu^ 1 in our final results. 

It should be pointed out that unlike the classical case, the pressure tensor does not necessarily vanish in the 
limit of zero temperature. However, we note that we need not be concerned about the contribution from the pressure 
term in this particular case. This is because Py vanishes linearly in the MHD limit and thereby enters as a cubic 
nonlinearity (which does not affect the three- wave interaction) in the evolution equation for Ejj. The pressure tensor, 
however, also gives a contribution in the MHD equation ©, but it turns out that this is a nonlinear contribution 
proportional to (u> cg — o; c ) _1 which is small compared to leading terms proportional to uj~ 1 . We may therefore neglect 
the contribution from P^- altogether. 
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One-fluid calculation 



We start by considering the one fluid model for which, as mentioned above, the unperturbed spin-density is zero, 
S = 0. Our first aim is to find the linear dispersion relation as well as the linear eigenvectors (polarizations) of 
the shear Alfven wave and the compressional Alfven wave. We note that in the MHD equation (J5]) we need an 
expression for the magnetization. We therefore begin by solving the spin-velocity evolution equation to find the 
S-tensor. Linearly, this is straightforward and we find 





k x B x uj c 


k z B y 




(«*-<) 
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k x ByUJ c 


k z B x 
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As can be seen, the components in (|12l) have different magnitudes, but we keep all of them at this stage in the 
calculation. Next we use the linear E-tensor in the spin evolution equation to find an expression for the linear 
spin S and thereby the linear magnetization M = fin^S. 
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(13) 



Here we have dropped contributions to components of M that are smaller by factors \ui 2 — u; 2 ff | /ui 2 g and/or u>/u> cg . 

Similarly to the classical 

(14) 



Substituting (fT3)) into (j6|) , the linear dispersion relations are obtained from ((6]) , ((7|) and 
ideal MHD case, the modes decouple into the shear Alfven wave described by 



£U(w,k) 



= uj 2 



k z C A 



1 - 



in 



w c«) 
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and the compressional Alfven (fast magnetosonic) wave with the dispersion relation 



D M s(w,k) 



k 2 



kyC A 



A . .. 



Note that the last terms in (fT^|) and ([15)1 are smaller than the first spin-modification in (TTS"]) as u 2 <C 



Furthermore, the linear eigenvector components for the shear Alfven wave are 
and 



0, Bi 
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v = — ~ v 



For the compressional Alfven wave (index MS) we instead obtain u y 

kf s B 



0, B% IS = 0, and 



B 



MS 



MS 



MS 



B 



[.MS R 

jus 

I.MS 

n - 1 - i, MS 
^ a ,,,MS * 



(15) 
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(16) 

(17) 
(18) 
(19) 



Next we aim to calculate the three wave coupling coefficients due to the quadratic nonlinearities. We have calculated 
the nonlinear contribution to the coupling coefficients including all terms proportional to {tu cg — wj" 1 . However, 
our results show that these terms only give rise to small corrections to the leading terms proportional to uj^ 1 . Since 
the full analysis is rather tedious we will therefore only write out the leading terms in the NL contribution to the 
S-tensor as well as to the magnetization M. Under the given approximations, keeping the resonant terms, we find 
the components with a nonzero nonlinear contribution to the S-tensor to be 



■lA' 
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(20) 
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and 



■lA' 



k x B z , i q k^ s B^Bf 



u r to cj r o; 



MS 



(21) 



for the Alfven wave. Here * denotes complex conjugation. For the magnetosonic wave the component with a nonzero 
nonlinear contribution is 



J 33 



k z B¥ s 2iu 1 
i 2 +-^ 

ui n Lo C gUj 



s {^ + kf)BfB* 



(22) 



Solving the spin evolution equation with the sources from E given above, we find the a nonlinear contribution to the 
^-component of the magnetization 



M MS = _ 



nofi 2 ( k\Bf s 2 M k f S ( k z + kA ' 



ft uj cg uj 2 ( MS ) v v 



(23) 



for the MS wave, and a nonlinear contribution to the y-component of the magnetization 
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(24) 



for the Alfven wave. Now that we have expressed the magnetization in terms of the magnetic field, correct to 
second order in the amplitude, we may substitute these results into ([5]), and perform the rest of the calculations 
using ([5]), (O and (JSJ as in standard MHD theory. Accounting for time dependent amplitudes with the substitions 
Da{uj,\l) — > [dDA/duj]id/dt and Dms{<jJ, k) — > [ODms /duj]id/dt, doing successive elimination keeping the velocity 
variables as the wave amplitudes, we find the following coupled equations for the different wave modes 



Ci, A \, MS 

dt "dD A ,/du y x 



(25) 



and 



with the coupling coefficient 



S)„.MS ,,,2(MS) 

-Cutut 



dt 



dD MS /d(j 



v y 



(26) 
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(27) 



Due to the symmetry between the two shear Alfven waves, the equation for dUy/dt is obtained by exchanging A 
and A' in Eq. (|2"S")) . The appearance of the common factor C in the three coupled equations is a reflection of the 
Manley-Rowe symmetry [25| ■ The first term of C is a purely classical contribution, that agrees with e.g. Refs. (28l [29| 
in the cold limit. For the spin contribution in Eq. (|27|) to be important as compared to the classical one, a rather 
dense plasma is required. By contrast, other MHD phenomena exists that require less extreme parameters for the 
electron spin to be important. Nevertheless, as will be discussed in the final section, the results derived here have a 
number of interesting theroretical consequences. It should be stressed that the contribution to the magnetization in 
this one-fluid model stems from the E-tensor. This is in contrast to the two-fluid model as we will see below. 



B. Two- fluid calculation 



We now consider the problem of three wave coupling using the two-fluid model. The spin is then determined from 
(|4"|) with the contribution from E omitted, as described in section II, but we now have two species of electrons, which 
have the unperturbed spin Sot = ^ an< ^ ^oi = ~ z, respectively. The total magnetization is then written as 



M = A(#S(t) +n (%W) 



(28) 
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which gives Mo = in agreement with the previous section, provided we let not = n oi = n o/2 which will be used 
henceforth. Next we find the linear spin-vector to be 



Si = 



-SoBy 







(29) 



Note here that although the terms oc So in (|2"9"|) are larger than the terms oc (uS^ — w^ g ) _1 , the former has opposite 
signs for the up- and down species, and hence give no contribution to the linear magnetization. It turns out that the 
terms in (|2"9"|) oc (o;^ — a;^ g ) _1 are needed to get agreement with the linear magnetization obtained with the one- fluid 
model (|13[) . However, it can be noted that these terms have been dropped in the expression for the coupling coefficient 
(|27p where, in the end, only the leading term is kept. Next we need to find an expression for the fluid densities of the 
electron spin fluids. From the continuity equation ([1} we have 

n« =n s) +n[ ) s) ^^+n^ (30) 

where ti( s ) nl is a non-linear contribution that can be shown not to contribute to the magnetization after summation of 
the spin states s ~ J [, An expression for the electron velocities is obtained by solving the fluid momentum equation 
@ together with — dtB = V x E. To close this system we may in general need to use a full fluid description. However, 
within the MHD regime and for this specific problem, it suffices to determine the magnetization. For our case with 
not = n oi = n o/2 several terms vanish in Eq. (|28p after the summation over up and down species. In MHD we make 
the approximation that we may write the electric field as E = — x B ~ u x B. This allows us to again make use 
of Eqs. (|16[) - (flT)|) . Under these assumptions we note that E z vanishes linearly, and also that E^ L is only proportional 
to quadratic combinations of B-field components and will therefore not contribute to the magnetization. Thus E z 
may therefore be set to zero from now on. Under these assumptions, it is easy to show that the fluid velocities may 
be written as 

v x = -. 1 [-v x B z + v z B x ) (31a) 

m uj c k x m oj c 

v y = - — — ^L + i-t!—k x B z S - — — (v y B z - v z B y ) - i-^—k x (B X S X - BySy) (31b) 



v z = — —k z B z S + i — — (v x B y - v y B x ) —k z (B X S X - B y S y ) (31c) 

mu! m to mui 



From Eqs. (|28p - (|3ip we find the linear magnetization, and it agrees identically with the expression obtained from 
the one-fluid model. Furthermore, an extended analysis gives agreement for the NL terms of the magnetization as 
well. Consequently, the coupling coefficients remain the same regardless if the one-fluid or two-fluid model are used 
to determine the magnetization. This corroborates the usefulness of the two fluid model in the MHD regime. 



IV. DISCUSSION 

In the present paper we have studied a recently presented fluid model accounting for the electron spin [16^ . and 
adopted it to the MHD regime. The main feature of the original model is that in addition to basic spin effects as 
the magnetic dipole force, spin precession, and the magnetization current it incorporates spin-velocity correleations. 
The spin-velcoity correlations have been shown to be important for a number of spin plasma phenomena (l6l . l2~i| . 
Introducing the approximations appropriate for the MHD regime, it turns out that essentially the ordinary MHD 
equations are recovered, but with a magnetization that needs to be determined. This can be done in different ways. 
Firstly from a single fluid model of electrons, that besides the basic spin precession contains spin- velocity correlations. 
Or, secondly, from a two-fluid model where spin-up and spin-down electrons constitute different species. A theoretical 
argument is presented in Section II B suggesting that these two models are equivalent in the MHD regime. However, 
the equivalence argument depends on certain assumptions which is difficult to justify rigorously, and thus practical 
tests of the equivalence is valuable. For this purpose we have evaluated the different models using a nonlinear three 
wave interaction as a test problem. Specifically, we have computed the coupling coefficients between two shear 
Alfven waves and a compressional Alfven wave. A classical as well as quantum mechanical (spin) contribution to the 
coupling coefficients are found, and the coupling coefficients are indeed identical in the two models. Furthermore, 
the coupling cofficients obey the Manley-Rowe symmetries [25[. The Manley-Rowe relations is a reflection of the 
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underlying Hamiltonian structure [3(3] of the model. The fact that the coefficients preserve these relations strongly 
suggests that the approximations made when deriving the models are sound, as otherwise it is highly likely that the 
Manly-Rowe symmetries would be broken. 

Since we have put forward two somewhat different models in this paper, one may ask which one that is most easy 
to use. For the analytical calculations made here, the degree of complexity is found to be roughly the same. However, 
the two-fluid model has a great advantage in case one would like to do Particle-In-Cell (PIC) simulations. In that 
case, the only modification of a standard code would be to have two species of electrons, and add a force proportional 
to ±/i e V-B in the momentum equation, as well as to compute M = |U e b(n-|- — rij,) to find the contribution from the 
Magnetization current in Amperes law. As the variables for the density and magnetic field are monitored throughout 
the PIC-simulations anyway, no new equations and only little extra complexity is added to the general concept. 
Developing a PIC-scheme incorporating spin- velocity correlations, however, is a much more cumbersome project, as 
the evolution equation for this object is not present in present PIC-schemes, and also such equations are considerably 
more complex. Furthermore, it is not clear that it is at all possible to model spin-velocity correlations as a single- 
particle property, which makes the adaption of this model to the PIC-scheme questionable conceptually. Thus we 
conclude that the two-fluid model is valuable for the purpose of incorporating electron spin effects in PIC-simulations, 
although we stress that the applicability of such an approach will be limited to the MHD regime. 
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